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ABSTRACT 

We explore the ability of observations of the global brightness temperature of the 21-cm 
signal to constrain the reionization history and the properties of the ionizing sources. 
In order to describe the reionization signal, we employ either a commonly-used toy 
model or a more realistic structure formation model that parameterizes the properties 
of the ionizing sources. If the structure formation model captures the actual evolution 
of the reionization signal, then detecting the signal is somewhat easier than it would 
be for the toy model; using the toy model in this case also leads to systematic errors 
in reconstructing the reionization history, though a sufficiently sensitive experiment 
should be able to distinguish between the two models. We show that under optimistic 
assumptions regarding systematic noise and foreground removal, one-year observations 
of the global 21-cm spectrum should be able to detect a wide range of realistic models 
and measure the main features of the reionization history while constraining the key 
properties of the ionizing sources. 
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1 INTRODUCTION 

One of the most important frontier fields of cosmology is the 
evolution of the Universe from the dark ages following hy- 
drogen recombination through to the epoch of reionization. 
The 21-cm line associated with the hyperfine transition of 
atomic hydrogen is the most promising signal for detecting 
and mapping the spatial and redshift distribution of hydro- 
gen in the universe, and for studying the sources responsible 
for heating and reionizing the intergalactic medium (IGM) 
at redshifts z J> 6. Indeed, an important feature of the 21- 
cm signal is that the spectral dimension allows in princi- 
ple 3D tomography of hydrogen as a function of redshift, 
providing much richer structure than the cosmic microwave 
background (CMB), which yields just a single sky map; this 
may he lp to detect primordial no n-Gaussianity and test in- 
flation |Loeb 8z Zaldarriagall20o3 ). During reionization, de- 
tecting the bubble structure would probe the main sources 
of ionizing radiation, even if these are otherwise unobserv- 
able because, for example, they are too faint to be detected 
individually. It is also important to characterize when and 
how long the reionization took place, because of the signif- 
icant effect of reionization on the subsequent formation of 
galaxies. In particular, the ultra-violet (UV) radiation heats 
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the gas, raising the Jeans mass and causing a suppression of 
star-forming galaxies in low-mass halos (< 1O 9 M ). 

A first estimate of the reionization redshift z r has 
been deduced from the CMB polarization, where an addi- 
tional peak on large angular scales, corresponding to the 
horizon size at the reionization epoch, is expected due to 
scattering, with an amplitude related to the total optical 
depth . The recent analysi s based on the 7-year WMAP 
data (jKomatsu et al.ll2010l ) finds a reionization redshift of 
10.4 ± 1.2, in terms of an equivalent instantaneous reion- 
ization, but is still consistent with a wide range of possible 
reionization histories. Other constraints come from the Lya 
galaxies at redshift 5.7 and 6.5, whose characteristic lumi- 
nosity function shows a lack of time ev olution that is consis- 
tent with a fully ionized IGM at z ~ 6 (jMalhotra fc Rhoads 
2004) . Analyses of the spectra of high-z QSOs (IFan et al 
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also suggest that the IGM is still very highly ionized at this 
redshift. 

The expected picture for reionization is thus an inho- 
mogeneous and extended process, for which the nature and 
the evolution of the ionizing sources are still observationally 
undetermined. Upcoming and future observational probes 
should allow us to distinguish among various reionization 
models and in particular constrain th e possible extended o r 
instantaneous nature of this process l|Bruscoli et al.ll2002l ). 
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Low-frequency observations with radio telescope arrays such 
as LOFAEQ, MWA0, PAPER0, and SKA0, will over the 
next decade constrain the spatial distribution of the ion- 
ization sources, while dipole observations of the spatially- 
integrated 21-cm signal are currently unde rway, although 
in their infancy, as exem p lified by EDGES {§ owman et al.l 
120081 ; iBowman fc RogerslhoiOTl 

Recently IPritchard &: Loebl l|2010l . hereafter PL) made 
a first theoretical attempt to predict detection limits for fu- 
ture observations of the global 21-cm signal. In this paper, 
we explore the potential for global 21-cm experiments to 
constrain reionization starting with the same simple ana- 
lytical toy model used by PL (based on the tanh function), 
but then focus on a more realistic, physically-based galaxy 
formation model that parameterizes the properties of the 
ionizing sources. In the following section we briefly review 
the general setup of global 21-cm measurements and the toy 
model, and then introduce a simple galaxy formation model 
within Cold Dark Matter (CDM) - dominated hierarchical 
structure formation. In section 3 we make predictions for the 
spatially-integrated 21-cm signal that experiments such as 
EDGES aim to measure. We compare the expected signal 
from the two models, and show some examples of the ex- 
pected errors in the model parameters that are reconstructed 
from observations. In section 3.4 we explore the systematic 
effects of assuming an incorrect model when trying to re- 
construct the global 21-cm signal. We finish this part with 
section 3.5, which presents our main result, the detection 
limits of the global 21-cm signal. We summarize and discuss 
our conclusions in section 4. 

Hereafter we assume a fiat A CDM cosmology, with 
matter density parameter Q. m = 0.272 (dark matter plus 
baryons), cosmological constant density parameter Qa = 
0.728, H = 70.4km/s/Mpc (Hubble constant), Q b = 0.045 
(baryons), ns = 0.963 (power spectrum index) and ug = 
0.809 (power spectrum nor malization) accordi ng to the lat- 
est 7-year WMAP results (jjarosik et al.ll2010f) . Unless oth- 
erwise stated, we estimate all errors at the 68.3 percent con- 
fidence level. 



2 MODELING THE 21-CM SIGNAL 

2.1 The 21-cm foreground and signal model 

In general, a global 21-cm measurement yields the antenna 
temperature T s k y (^) = Tf s (u) + T b (v), where Tf g (v) and 
T b (u) are the foreground and cosmological 21-cm brightness 
temperatures, respectively. For the cosmic signal, we assume 
that the dipole antenna temperature essentially measures a 
sky average, since fluctuations are expected to be present 
only on angular scales that correspond to small fractions of 
the sky. The foregrounds (i.e., our Galactic emission and ra- 
dio emission from other galaxies) have large-scale angular 
structure, but even if they are convolved with an angular 
dipole response, this does not affect our analysis, which only 
assumes that they are smooth as a function of frequency. 
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For the foreground brightness temperature Tf g , we as- 
sume a polynomial fit of the form 



log Tfg = a i lo &( v l 



(1) 



In particular, we use the third order polynomial fit from 
PL, who fitted the m o del o f the sky put together by 
Ide Oliveira-Costa et al.l (|2008T l using all existing observa- 
tions, by averaging the foregrounds over the dipole's angular 
response: 



log Tf g = log T + ai log(^/^o)+a2[log(^/i/o)] 2 

+ a 3 [log(^/^ )] 3 



(2) 



with parameter values vq — 150 MHz, To = 320 K, ai = 
—2.54, d2 = —0.074, and 03 = 0.013, chosen from fitting 
to the band v = 100 - 200 MHz. Note that at these fre- 
quencies Tfg is dominated by diffuse synchrotron radiation 
from the Galaxy. The residuals related to such a parame- 
terization of the foregro und are dominated by limitatio ns of 
the adapted sky model |de Oliveira-Costa et al.ll2008l ) and 
they are < 1 mK averaged over the band. While in princi- 
ple higher order polynomials may be needed to reduce such 
residuals in the future, given the smoothness of the spec- 
trum of the foreground, low order polynomials are key to 
avoid throwing the signal away with the foreground and to 
reduce the statistical errors (§ 13. 5[) . 

For the cosmic 21-cm signal, the brightness temperature 
through the IGM is T b = T C MBe" T + T s (l - e~ T ), with 
t(z) <C 1 the optical depth at 21(1 + z) cm produced by a 
patch of neutral hydrogen at the mean density and with a 
uniform 21-cm spin temperature Ts, 



r(z) = 9.0 x 10~ 



TcMB \ / flft/A ( flm 



T s A0.03A0.3 



10 



(3) 



During the epoch of reionization the Lyman-a and X-ray 
radiation backgrounds are expected to be strong enough 
to bring the spin temperature Ts to the gas temperature 
and heat the cosmic gas well above the cosm ic microwave 
background temperature |Madau et al.|[l997l ). Under these 
conditions, the observed 21-cm brightness temperature T b 
relative to the CMB temperature Tcmb is independent of 
Ts. Therefore, T b (hereafter measured relative to Tcmb) is 
given by: 



T b {z) =(T s -TcMB)(l-e- T )Q H i = 
1/2 

=T 21 (^ T7 ^) Qui 



10 



(4) 



where T21 = 9.0 x 10- 3 (fi(,/z/0-03)(n m /0.3)- 1/2 T C mb 
27.2 mK, and Qui = Ahi/(A H i + Ami) is the neutral hy- 
drogen fraction. Note that the ionized fraction is Qu 11 = 
1 — Qhi- Throughout this paper, given that we are inter- 
ested in the spatially-integrated 21-cm signal, we consider 
only the cosmic mean neutral or ionized fraction, and ne- 
glect spatial fluctuations in the 21-cm signal from density 
and peculiar velocity fluctuations. 

We consider an experiment covering the frequency range 
100 - 250 MHz in 50 bins of bandwidth B = 3 MHz for 
each of the receiver frequency channels, and integrating time 
tint = 500 hours (these parameters mimic EDGES with an 
order of magnitude longer integration time). Under these 
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assumptions, the thermal noise in the i'th receiver frequency 
channel is given by the radiometer equation: 



2 T Bk (t>i) 

<7, = 



Btm 



(5) 



We note that the frequency range we consider corresponds 
to the redshift range 4.7-13.2. 

Our model thus consists of the foreground brightness 
temperature Tf g (y) and a suitable model for the cosmolog- 
ical 21-cm signal T{,{y). To derive the parameter errors, we 
directly calculate the Fisher matrix of the foreground plus 
21-cm signal parameters p expected with the above thermal 
noise ai, 



E 



J_ / <9r sky (u n ;p) dT sky (u„;p) 



dpi 



dpj 



(6) 



This equation provides an estimate of the covariance matrix 
C = F -1 , and therefore of the parameter uncertainty in 
dipole observations. Note that this is equivalent to finding 
the covariance matrix near the minimum \ 2 ■ These errors 
should be accurate as long as they are small. However, in 
many cases we consider regions of parameter space where 
the errors are large, e.g., when we calculate the detection 
limit of an experiment, or more generally due to parameter 
degeneracies. Thus, we often use a more generally-valid and 
computationally intensive Monte-Carlo (MC) error analysis. 
We generate a large number of MC simulations of the mea- 
surement noise, finding the best-fit parameters in each case 
by minimizing the \ 2 '- 



E 



[AT 3ky (u n ;p)Y 



(7) 



where AT s k y is the difference between the measured and pre- 
dicted total 21-cm sky temperature in channel n (centered at 
the frequency v n ). The distribution of best-fit parameters in 
the MC trials yields parameter errors and their correlations. 

In the following sections we will focus on the modeling 
of the 21-cm signal, in particular of the neutral fraction Qu i 
in equation Q. Given the great uncertainty associated the 
evolution of the neutral fraction Qh i due to the uncertain 
astrophysics of the ionizing sources, we will begin with a 
toy model, namely the tanh-based parameterization used by 
previous authors, which simply characterizes when and for 
how long the reionization occurs; then we will consider a 
more complex and physically motivated structure formation 
model, in order to better describe the reionization process 
and extract interesting astrophysical information, such as 
the mass of the smallest galaxies that can form and con- 
tribute to the IGM ionization, the overall number of ioniza- 
tions per baryon and the redshift evolution of the ionizing 



2.2 The tanfo-hased model of reionization 

The tanh-based parameterization is characterized by two pa- 
rameters describing the two main features of reionization: its 
mid-point z r and duration Az. This approach was used b; 
PL for the 21-cm signal (note that I Bowman et al.l (1200. 
used a somewhat different parameterization), and a similar 
tanh-based fitting function is the default parameterization 
of reionization in CAMB (although there it is based on the 



optical depth for CMB scattering) (|Lewisll200sh . Under the 
assumptions outlined above for the gas state during reion- 
ization, the 21-cm signal is given by 

T b (z) = T 21 (l±f\ \ [tanh (*^) + ij . (8) 



10 



Note that z r is the redshift at which the ionized fraction 
Qu ii = 50%, while z r + Az and z r — Az are the redshifts at 
which Qhii = 11.9% and 88.1%, respectively. This param- 
eterization is a convenient mathematical toy model but it 
does not have any particular physical motivation. We con- 
sider both the case where we fix the amplitude of the signal 
T21 to its known value (equation [4}, and the case where we 
leave it as a free parameter (following PL). 



2.3 A simple CDM-dominated galaxy formation 
model 

In the previous section we considered a toy model that has 
been used in previous observational and theoretical papers. 
While a toy model can be justified as an unbiased anal- 
ysis tool, especially given the large current uncertainty in 
the astrophysics of high-redshift galaxies, such an approach 
is also problematic. The particular model assumed (with a 
fixed, arbitrarily-chosen shape) may lead to systematically 
biased results if it cannot reasonably approximate the real 
reionization (we consider this issue further below). In addi- 
tion, it can be hard to interpret any results of a toy model 
in terms of the underlying parameters of interest. In partic- 
ular, the redshift evolution of reionization is closely related 
to structure formation. Indeed, reionization is driven by the 
intergalactic ionizing radiation field, which (we expect) is 
the result of the ionizing radiation escaping from stars and 
quasars within galaxies. While astrophysical aspects (such 
as star formation and feedback) play a significant role, the 
evolution of galaxies is driven by the properties of the host 
dark matter halos. A major reason for studying reionization 
is to learn more about both galaxy formation and the astro- 
physical properties of galaxies in the reionization era. Thus, 
a more realistic and useful approach is to use models based 
on our understanding of CDM-driven galaxy and structure 
formation, a model with many successes at lower redshifts, 
and to include some flexibility in order to account for the 
uncertain astrophysical parameters. Here we take the first 
step in this process by using a simple model that is based 

on the standard theory of galaxy formation. 

W e begin with the equation from iBarkana fc Loebl 
(2001), based on lShapiro fc Girouxl l|l987h . that statistically 
describes the transition from a neutral universe to a fully 
ionized one; in particular this equation describes the evolu- 
tion of the HII filling factor Qhii, i.e., the fraction of the 
volume of the universe which is filled by H II regions. 



dQ H 11 JVion dF a 



C -0 r\ 
aB^rn H Qu 11 



dt 0.76 dt ~" a 3-- ' ( 9 ) 
assuming a primordial mass fraction of hydrogen of 0.76. In 
this equation Ni on = iV 7 / s t ar /esc is an efficiency parameter 
that gives the overall number of ionizing photons per baryon; 
for instance, if we assume that baryons are incorporated into 
stars with an efficiency of /star = 10%, the escape fraction for 
the resulting ionizing radiation is / osc = 5% and N y m 4000 
ionizing photons are produced per baryon in stars (for a 
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stellar IMF similar to the one measured locally but with a 
metallicity equal to 1/20 of the solar value), we infer that for 
every baryon in galaxies ~20 escaping ionizing photons are 
produced by stars. We obtain a similar result if we co nsider 
mini-quasars rather than stars (|Barkana fc Loeb| [2001). It is 
possible to get a substantially higher 7Vi on using Pop III stars 
or by assuming a high escape fraction. Ni on also determines 
the maximum comoving radius of the region that a halo of 
mass M can ionize on its own (neglecting recombinations), 

W = 675^0^^) , (10) 

a radius that is larger than the halo virial radius by a factor 
of ~ 20 (essentially independent of redshift and halo mass). 

Also in equation (J5J), a = 1/(1 + z) is the scale factor, 
n% is the present number density of hydrogen, as is the 
case B recombination coefficient of hydrogen, and C repre- 
sents the volume-averaged clumping factor (in general time- 
dependent), 

C=(n 2 H )/n 2 H . (11) 

This factor crudely accounts for a non-uniform IGM that 
includes high-density clumps. Since each ionized bubble is 
far larger than the typical scale of clumping, so that many 
clumps are averaged over, C can be assumed to be approx- 
imately spatially uniform. 

The collapsed fraction F co \ is the fraction of all the 
baryons in the universe that is in galaxies, i.e., the fraction of 
gas which settles into halos and cools efficiently inside them. 
A simple estimate of the collapse fraction at high redshift is 
the halo mass fraction above some cooling threshold. More 
generally, we include halos above some minimum circular 
velocity V c . We use the Sheth-Tormen halo mass function, 
which accurately fits the me an halo abundance in simula- 
tions (|Sheth fc Tormenll2002l '). We calculate the power spec- 
trum transfer function using the CAMB code fL ewis et al.l 
|2000D . 

The solution of equation is (Barkana & Locb 2001) 

r\ u\ f t Ni on dF co \ F(t't),,l , 1n \ 

QHIlW = y OT6^ e dt - (i2) 

where (if C is time-independent) 

F{t',t) = -\^0L-C[f{t')-f{t)\ , (13) 
and where (in flat ACDM) 

Once Qhii(£) reaches unity, the universe becomes fully 
reionized and remains so within our model. 

Equation (|12[) allows us to quickly calculate the time 
evolution of the ionized fraction of the universe once we fix 
the IGM clumping factor and the parameters related to the 
ionizing sources. Also, in equation ((6| we calculate accurate 
derivatives for the 21-cm signal as numerical integrals of par- 
tial derivatives of the integrand in equation Hereafter, 
we refer to this CDM-dominated galaxy formation model as 
the CDM model. 

Within the CDM model, the parameters that determine 
the redshift evolution of Qhii (and Tb) are JVi on (which 



we assume is a constant, in this first investigation of fit- 
ting global 21-cm signals from a galaxy formation model), 
C (likewise assumed constant), and the minimum halo cir- 
cular velocity V c (equivalent to a minimum mass) required 
for halos that host galaxies. We allow V c to vary, since while 
cooling sets a minimum value for it, feedback (radiative or 
from supernovae) may in reality set a higher threshold for 
effective star formation. We set C — 1 as our standard value 
(i.e., corresponding to a uniform IGM), and discuss in sev- 
eral cases the effect of allowing C to vary. 

3 RESULTS 

3.1 Global properties of the two models 

We begin by visually comparing our two models, the tanh- 
based toy model and the more realistic CDM-based model. 
In Figure [T] we plot a few examples of the global (volume- 
averaged) 21-cm signal for each model, over the assumed ex- 
perimental frequency range (100-250 MHz). For the CDM 
model, we consider a minimum halo circular velocity V c = 
{4.5,16.5,36,64} km/s, corresponding to a midpoint of 
reionization z r — {11.8,8.7,6.9,5.6} and a minimum galac- 
tic halo of 1.3 x 10 6 , 1.0 x 10 s , 1.4 x 10 9 , 1.0 x 10 10 M at 
z — z r , with the other parameters fixed at iVi on = 20 and 
C — 1. The two lower values of V c correspond to cooling 
via molecular hydrogen and via atomic hydrogen and he- 
lium, respectively, while the higher values are roughly in the 
range of values possible due to photoheating or supernova 
feedback. In general, a higher circular velocity implies that 
only more massive halos are included, delaying reionization 
to a lower redshift (for a fixed Ni on ). In addition, since more 
massive halos are more rare and correspond to the Gaussian 
tail of large (positive) density fluctuations, their abundance 
changes rapidly with time, thus making reionization more 
rapid in terms of its redshift extent; this would be obvious in 
a comparison of models all normalized to the same z r , while 
the effect is suppressed in this figure since the lower redshifts 
are stretched into a relatively large frequency interval. The 
main point is that z r is determined by a combination of V c 
and A r i on , while the extent of reionization is separately sensi- 
tive to V c , so that a measurement of the reionization history 
can probe both the characteristic halo mass of the ionizing 
sources and their ionizing efficiency. 

The main qualitative difference between the two models 
is that the CDM model shows a steady rise of Qhii, while 
the toy model is much more round in shape, in particular 
showing a slowdown of reionization during its last quarter or 
so. The toy model is explicitly symmetric in redshift about 
the midpoint of reionization, while in the CDM model reion- 
ization starts slowly but ends quickly. The steady acceler- 
ation of reionization in the CDM model is driven by the 
exponential rise of the ionizing sources, which correspond to 
rare halos at these redshifts. However, our simple model is 
by no means fully general, so we treat our conclusions with 
caution, as discussed further below. We consider the CDM 
model to be an example of a realistic model, which may be 
quantitatively plausible if some of the missing complications 
turn out to have a relatively minor effect on the global 21- 
cm signal. We note, however, that some complications may 
tend to make the CDM model more similar in shape to the 
tanh model (see Sect. 4). 
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Figure 1. Global 21-cm signal predicted by the CDM model 
(solid curves) for various values of the minimum circular velocity 
of galactic halos (V c = {4.5,16.5,36,64} from left to right). We 
fix C = 1 and N[ on = 20. We also show (dashed curves) the 
global 21-cm signal for a tanh model of reionization with Az = 2 
and z r = {6,8,10} (from right to left). The signal for a fully 
neutral universe is shown for comparison (dotted curve). Note 
some redshift values (at the top) that fall within the experimental 
frequency range (z = 4.7 — 13.2). 



The tanh-based model is explicitly expressed in terms 
of the mid-point z r and duration Az of reionization, while 
in the CDM model these are derived parameters. While the 
midpoint z r is naturally defined as Qu n = 50%, there is 
some ambiguity in Az. For the toy model, we have chosen 
to follow previous PL in defining Az as above, a definition 
that is natural for the tanh function, and implies that z r + 
Az and z r — Az (i.e., a total spread of 2Az) delineate the 
central 76.2% of reionization. However, for the CDM model 
we use a definition that should be the natural one more 
generally: Az — (z_io- — z + i CT )/2, with z_i CT and z + i a being 
the redshifts corresponding to Qhii = 0.16 and Qhii = 
0.84, respectively. Thus, in the CDM model a spread of 2Az 
marks the central 68%. In the context of the tanh model, 
this definition would give a value of Az smaller by a factor 
of 1.2 than the definition that we have followed. 

In order to gain intuition on how the characteristics 
of reionization are set in the CDM model, in Figure [5] we 
show the dependence of z r and Az on JVi on for one value 
of V c and several values of C . Larger values of JVi on lead to 
earlier reionization (i.e., higher z r ) at a time when the ion- 
izing sources are brighter and rarer, so their rarity leads to 
a shorter span Az for reionization. We compare C = 1 to no 
recombinations C — and fast recombinations (C = 10). At 
least during most of reionization C is likely to be of order 
unity, since the low-density IGM gets reionized first, and the 
denser gas is left for the final stages of reionization. We find 
that a high clumping factor can be essentially counterbal- 
anced by a higher value of iVi on , at least during the central 



Figure 2. Reionization characteristics z r and Az as a function 
of Ni on for the CDM model. We fix V c = 16.5 km/s and consider 
C = {0, 1, 10} (dotted, solid and dashed curves, respectively). 



portion of reionization that defines z r and Az, so that in- 
cluding C as a free parameter mostly includes the degener- 
acy of the parameters but does not significantly change the 
allowed parameter space of z r and Az. 

A more complete picture of the allowed parameter space 
is shown in Figure O where we present the isocontours of z r 
and Az in the Ni on -V c plane. For reasonable values of AW 
and Vc the reionization midpoint z r varies widely (from be- 
low 6 to above 18), while the span of reionization Az covers 
roughly the range 1-3. As noted above, these values of Az 
should be multiplied by a factor of 1.2 for a fair comparison 
with Az in the tanh model. 

3.2 Expected parameter errors: the tanh model 

We now derive parameter errors for some specific instances 
of potential global 21-cm observations. We begin with the 
tanh model, and test the Fisher matrix formalism against the 
MC error analysis. Taking fiducial values of z r — 8, Az — 1, 
and assuming the model parameters for the foreground as in 
§ 12.11 we generate 10 MC simulations of the noise, finding 
the best fit parameters in each case. The resulting parameter 
contours are shown in Figure [3] along with the corresponding 
Fisher matrix constraints. We find good agreement between 
the two methods in this example. 

We consider both the case where T21 is a free parame- 
ter (as assumed by PL) and where it is fixed at its known 
value. The error ellipses show that there is a strong positive 
correlation between T21 and Az, i.e., there is an uncertainty 
in distinguishing between a higher amplitude extended sce- 
nario and a lower amplitude quicker scenario, since both 
produce a similar slope with frequency in the 21-cm signal, 
and it is this sharp slope that can be distinguished from 
the foregrounds (which are smooth and thus can be mod- 
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Figure 3. Isocontours of z r (solid lines) and Az (dashed lines) 
derived from the CDM model in the N^ on -V c plane (for C = 
1). The parameters (Ni on ,V c ) are allowed to vary over iVj on £ 
{2, 1000} and V c g {4.5, 100} km/s. 
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Figure 4. The 68 and 95% confidence regions of various pa- 
rameter pairs for the tank model of reionization, comparing the 
MC likelihood (green/bright and red/dark shaded regions, respec- 
tively) to the Fisher matrix (solid ellipses) calculations. For com- 
parison with PL, we set z r = 8 and Az = 1, T21 free and fit four 
polynomial (foreground) parameters (i.e., N po i y = 3). In the pan- 
els on the right, we also plot the Fisher matrix (dashed ellipses) 
results for the same model with T21 fixed at its known value. We 
assume an integration time tj nt = 500 hours. 



eled by a low-order polynomial). There are also significant 
correlations among the other parameters. Comparing with 
PL, we note that the amplitude of the errors that we find are 
smaller in all four panels, and also the sign of the correlation 
is different in the z r — T21 and z r — Az relations. Note that 
PL used a different frequency interval, i.e. 100 — 200 MHz 
(Pritchard & Loeb, personal communication). 

The results thus show significant correlations among the 
parameters when T21 is free, but substantially reduced er- 
rors and correlations when T21 is fixed. We conclude that it 
is possible to obtain a direct observational estimate for T21 
from these type of data, in order to check consistency with 
the theoretically expected value, but in order to constrain 
reionization it is very helpful to use our independent knowl- 
edge of T21. In this example we have assumed an integration 
time tint = 500 and fitted a foreground polynomial of degree 
A^poiy = 3 over the entire frequency range, assuming no re- 
maining foreground or systematic residuals. This represents 
a quite optimistic assumption regarding the level of system- 
atic noise and the ease of foreground removal, far beyond 
the current EDGES experiment, as discussed further below. 

Since we have just considered a rather optimistic exper- 
imental scenario, it is interesting to consider more realistic 
possibilities. One way to do this is to vary the integration 
time, thus increasing the errors. We can take this also as a 
rough indication of the effect of increasing the foreground or 
systematic residuals to various levels (still with iVpoiy = 3). 
In the case considered in Figure [H the errors per frequency 
bin range from 0.4 mK in the lowest-frequency bin to an or- 
der of magnitude lower at the highest-frequency bin. More 
generally, the noise varies with the integration time tint of 
the bolometer as t^' 5 [equation J5J], In Figure [5] we con- 
sider the fractional error on z r as a function of tint for 
Az = {1,2,3} and z r = 8. Note that the fractional error 
in z r varies approximately as tjnt 5 since the errors are rel- 
atively small over most of the plotted range (which makes 
the model behave approximately like a linear model). More 
extended reionization scenarios increase the errors signifi- 
cantly. Fixing T21 at its known value reduces the errors by 
15 - 50%. 

We will directly consider detection limits in a later sec- 
tion, but one way to define a successful detection of reioniza- 
tion is when observations yield a meaningful constraint on 
the most interesting single number associated with reioniza- 
tion, namely z T . Within the tanh model, rough (10%) con- 
straints on z r are expected for tint = 26 hours (if Az — 3) or 
1.9 hours (if Az — 2), while tight (1%) constraints require 
tint = 848 hours (if Az = 3), 51 hours (if Az = 2), or 3.1 
hours (for sharper reionization, with Az — 1). 

In Figure [6] we show a different range of the parameter 
space, considering three possible values of z r while varying 
Az over a wide range, all for tj n t = 500 hours. Here we show 
the relative errors on both Az and z r , finding that z r is gen- 
erally better constrained, by up to an order of magnitude. 
The errors increase with Az, roughly saturating at 30% for 
z r and 50% for Az (i.e., the errors only increase slowly be- 
yond these values as Az is further increased beyond ~ 5). 
As before, fixing T21 at its known value can make a big 
difference (compared to allowing it to be a free parameter) , 
especially in constraining Az (except when all the errors are 
large, for high Az). The relative errors vary weakly with z r 
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t (hours) 

Figure 5. Relative error on z r as a function of ti nt for Az = 
{1, 2, 3} (from bottom to top in each case) and z r = 8 for the 
tanh model of reionization. We consider T21 fixed or free in the 
analysis (solid and dashed curves, respectively). The errors have 
been calculated via MC analysis. 



over the range of 6-10 (a range which is all well within our 
assumed experimental frequency window). 

3.3 Expected parameter errors: the CDM model 

We begin our exploration of global 21-cm measurements in 
the context of the CDM model with Figure[7] where we show 
parameter errors and correlations for the fiducial values of 
(20, 1, 16.5) for (AW, C, V c ), respectively. The error ellipses 
show that there is a strong positive correlation between V c 
and iVi on ; indeed, this is a partial degeneracy, since while the 
error ellipse covers a small total area, each of these param- 
eters is uncertain at a relatively high (~ 10%) level. From 
Figure [3] it is apparent that this degeneracy with a positive 
correlation is driven by the value of z r , which is the main 
constraint from these observations (at least in the example 
we are considering of a high-precision experiment with low 
noise) . There is also a strong anti-correlation between V c and 
To — (To), demonstrating how the foreground fitting removes 
power from the total signal T s k y (y) making it more difficult 
to determine the parameters of reionization. The Fisher ma- 
trix and MC formalisms yield a reasonable agreement, but 
there are bigger differences compared to the tanh model, 
likely because the partial degeneracy gives larger errors in 
some directions in the CDM model. 

Fortunately, the partial degeneracy in the parameters of 
the CDM model is relatively harmless in terms of measuring 
the characteristics of reionization. For the case considered 
in Figure [7] we measure z r = 8.74 ± 0.02 and Az = 1.83 ± 
0.02. The two-parameter contour is shown in Figure [H] The 
relative errors in z r and Az are much smaller than in V c and 
Mon, showing that the global 21-cm measurements constrain 




12 3 4 5 6 7 8 
Az 

Figure 6. Relative error on Az and function of Az for 

the tanh model of reionization, for z r = {6, 8, 10} (dotted, dashed 
and solid curves, respectively). In each case we consider T21 to 
be fixed or free (where fixed corresponds to the lower curve at 
the left end of the plot). The errors have been calculated via MC 
analysis. 
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Figure 7. The 68 and 95% confidence regions of various parame- 
ter pairs for the CDM reionization model, comparing the MC like- 
lihood (green/bright and red/dark shaded regions, respectively) 
to the Fisher matrix (solid ellipses) calculations. We set our model 
parameters (JVj on , C, V c ) to the fiducial values of (20, 1, 16.5) and 
assume an integration time t- lnt = 500 hours and a foreground 
(plus systematics) polynomial with N po i y = 3. 
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Az 

Figure 8. The 68 and 95% confidence region of the reioniza- 
tion characteristics z r and Az derived from the CDM model, 
comparing the MC likelihood (green/bright and red/dark shaded 
regions, respectively) to the Fisher matrix (solid ellipses) calcu- 
lations. We set our model parameters (N[ on , C, V c ) to the fidu- 
cial values of (20, 1, 16.5) and assume t- lnt = 500 hours and 
A^poly = 3. The black dot corresponds to the input model val- 
ues («r, A*) = (8.74, 1.83). 



these quantities rather directly, somewhat independently of 
the underlying galaxy and halo parameters. Also, as noted 
above, the fractional error on z r is significantly smaller than 
in Az. The plotted results assume (7 = 1, but we find that 
if we allow the clumping factor to be a free parameter in the 
fit, this increases the CDM model parameter degeneracies 
but it does not significantly affect the errors on z r and Az. 

As we did for the tanh model, we vary the integration 
time and consider the expected experimental accuracy in 
measuring the most important quantity, z r . Figure [5] shows 
the fractional error on z r for Az — {1.5,2,2.5}. As in Fig- 
ure[5l the fractional error varies approximately as t^' 5 , and 
increases for larger Az. For C = 1, tight (1%) constraints 
on z r require tint = 68 hours (if Az — 2.5), 29 hours (if 
Az = 2), or 13 hours (if Az = 1.5). The CDM model gives 
somewhat better accuracy than the tanh model, for similar 
values of z r and Az, though the numbers are comparable. 



3.4 Systematic effect of the choice of reionization 
model 

Our use of two different models allows us to explore the sys- 
tematic effects of assuming an incorrect model when trying 
to reconstruct the global 21-cm signal from observations. 
We assume our more realistic CDM model as the input 
model, and try to fit the resulting 21-cm signal with the 
tanh model. In Figurc[]IJ]we plot the 21-cm signal as inferred 
from the fit of the tanh model + foreground to the 21-cm 
data generated from the CDM model + foreground with 



F 




1 10 100 1000 10000 

t (hours) 

Figure 9. Relative error on z r as a function of t- lnt for a CDM 
model of reionization with Az = {1.5, 2, 2.5} (black, red and blue 
colors, respectively; also from bottom to top in each case). We set 
z r = 8 and consider C = 0, C = 1 and C = 10 (dotted, solid and 
dashed curves, respectively). The errors have been calculated via 
MC analysis, with C held fixed in the fitting. 



(Ni on , C, V c ) = (20, 1, 16.5), corresponding to z r — 8.74 and 
Az = 1.83. The fit of the tanh-based model + foregrounds, 
after the subtraction of the best-fit foreground polynomial 
(which takes out part of the signal together with the original 
input foreground), leads to a quite different output profile 
of Tb(v) compared to the input one, and to biased values 
of the midpoint and duration of reionization. The best-fit 
parameters are z r = 8.19 ± 0.01 and Az — 1.28 ± 0.01; the 
latter is even more discrepant than may appear, because the 
input CDM value of Az should be multiplied by 1.2 for a 
fair comparison with the tanh model. While the statistical 
errors of the fit are tiny (for ti nt = 500 hours and N po \ y = 3), 
the systematic errors are quite large. The systematic errors 
are related both to the inadequacy of the tanh model in rep- 
resenting the reionization signal and to the presence of the 
foreground which must be fit with a polynomial. 

The news, though, is not all bad, since an experiment 
with such low noise levels would result in high, strongly dis- 
crepant, x 2 values for such a poor fit, giving a clear indica- 
tion that the template being used must, indeed, be modified. 
In particular, we find \ 2 ~ 1044 for 44 degrees of freedom. 
This means that in this example, only a much reduced exper- 
imental sensitivity corresponding to tint ~ 20 hours would 
give a reduced \ 2 °f order unity for the tanh model fit. Any 
experiment above this sensitivity would be able to discrim- 
inate between the CDM and tanh models. 

In Figure [11] we explore these kind of systematic errors 
over a wider range of the parameter space. We compare the 
best-fit parameters z r and Az from fitting the tanh model 
+ foreground with the true values for an input CDM model 
of reionization. We fix the input Az = 1.5 (equivalent to 
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Figure 10. Fit of a tanh model + foreground to 21-cm data gen- 
erated from a CDM model + foreground. We show just the 21-cm 
signal (without the foreground), for the input CDM model (points 
with error bars showing the measurement errors) and for the best- 
fit tanh model (solid curve). The input model sets (Ni OIL , C, V c ) 
to the fiducial values of (20, 1, 16.5). We obtain z r = 8.19 ± 0.01 
and Az = 1.28 ± 0.01 as best fit parameters with \ 2 = 1044 (for 
44 degrees of freedom), while the true values for the CDM model 
are z r = 8.74 and Az = 1.83. We assume t- lrLt = 500 hours and 
Npoly = 3. 



~ 1.8 in the tanh model), and vary z r over the range 6-10. 
This wider range shows similar results to the example above, 
where the best-fit z r in the tanh model is underestimated 
typically by 5-10%, while Az is underestimated by much 
more. We consider several different values of the assumed 
clumping factor in the CDM model, and find that the C = 
and C — 1 curves lie on top of each other, while C — 10 is 
only slightly different (where the comparison is made with 
fixed z r and Az values in the input CDM model). 



3.5 Detection limits of the global 21-cm signal 

In this section we present our main result, i.e., the experi- 
mental sensitivity that is required to detect the global 21- 
cm signal, as predicted by each of the reionization mod- 
els. A range of different results is summarized in Figure [T2l 
First we display the full range of allowed values of the mid- 
point and span of reionization within the CDM model (gray 
shaded region), where the parameters (Ni on , V c ) are allowed 
to vary over AT ion g {2, 1000} and V c G {4.5, 100} km/s, 
fixing C = 1. This region reflects Figure [3] showing that a 
wide range of z r is plausible, while the most relevant range 
of z r = 6 — 12 includes some models with Az as low as ~ 1. 

The figure also shows curves which delineate the 95% 
detection region for the different reionization models, for 
various polynomial orders of the fit (iV po i y = {3, 6, 9, 12} in 
equation [T]), and several possible values of the integration 




6 7 8 9 10 



Figure 11. Values of z r and Az from the fit of a tanh model 
+ foreground to 21-cm data generated from a CDM model of 
reionization with Az = 1.5, as a function of the input z r of the 
CDM model. We consider C = 0, C = 1 and C = 10 (dotted, 
solid and short-dashed curves, respectively); the C = and C = 1 
cases are indistinguishable. We also indicate the input values for 
the CDM model (long-dashed curves). T21 has been kept fixed 
to its known value in the fits. We have assumed AT po i y = 3. This 
plot is independent of tj nt . 



time ti nt . We define a model as detected if it is inconsistent 
with a fit that does not include the reionization signal. Thus, 
we fit each model signal with a foreground polynomial of a 
particular N po \ y , and if the resulting minimum x 2 is incon- 
sistent with zero at greater than 95% confidence, than that 
model is included within the detection region. For the CDM 
model, the unrealistic sudden end to reionization raises the 
X 2 values somewhat, so we reduce our sensitivity to this by 
removing from the y 2 value the contribution of the 12 MHz 
band centered on the end of reionization for each model. 
For both models there are some oscillations in the detection 
limit curves due to the degeneracy between the 21-cm signal 
and the polynomial fitting of the foreground. 

For the tanh model, for the same case as PL (ti nt = 500 
hours) we find significantly better prospects for detectabil- 
ity, with our limits on Az for a given z r typically higher 
by a factor of ~ 1.5 than their result. It is possible to con- 
strain some models with Az ~ 1 even with t- lnt = 10 hours, 
if iVpoiy = 3 suffices for removing the foreground and other 
systematics, or with tint = 500 hours if N po i y = 9. The worst 
case of iVpoiy = 12 only allows the detection (or ruling out) 
of very sharp reionization models that are probably unreal- 
istic. 

Our CDM model gives comparable constraints to the 
tanh model for low values of -/V po i y , but it is significantly 
more detectable with higher-order polynomials. The rapid 
rise of reionization up until its sharp end is easier to dis- 
tinguish from a high-order polynomial compared with the 
smooth tanh model, even after the removal (in the CDM 
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Figure 12. The 95% detection region for global 21-cm experiments, in terms of the midpoint z r and span Az of reionization. We consider 
polynomial order in the fit N po i y = {3, 6, 9, 12} in the various panels. We consider both the CDM model (solid curves) and the tanh 
model (dashed curves), in each case for an observational integration time ti nt = 500,50, or 10 hours (red, blue, and green, respectively, 
also from top to bottom). We also show the full range of allowed values of z r and Az of reionization (gray shaded area) assuming the 
parameter space 7V ion e {2, 1000}, V c e {4.5, 100} km/s, (7 = 1. 



case) of the frequency interval right near the end of reion- 
ization. We thus find that an interesting parameter space 
of CDM models can be detected even with iVpoiy = 12, for 
integration times of at least ~ 50 hours. 

We have not tried to indicate current constraints on 
reionization in in Figure [T2l to avoid overcrowding the figure, 
especially since these constraints are not directly expressed 
in terms of z r and Az, and the conversion to these variables 
would differ somewhat between the tanh and CDM models. 



Roughly, for these models, the 7-year WMAP data implies 
a 95% limit of z r > 8, while the absorption constraints that 
show a high ionization fraction at z ~ 6.5 imply a minimum 
z r that increases beyond 8 if Az > 1 (see also the Intro- 
duction and the discussion in PL). The relation between 
these constraints on reionization and those from global 21- 
cm measurements would change for more complex models of 
reionization. 
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4 CONCLUSIONS 

The aim of this paper was to investigate the possibility that 
global 21-cm observations during the epoch of reionization 
can probe the evolution of the IGM and the physical proper- 
ties of the ionizing sources. Detecting the 21-cm signal in the 
presence of the large foregrounds is challenging and it is im- 
portant to explore all avenues. While interferometric radio 
arrays are gearing up to measure 21-cm fluctuations, global 
measurements with a single dipole experiment can provide 
an independent and complementary method for detecting 
and/or constraining reionization. 

In order to derive quantitative predictions, we have im- 
plemented both a previously- used toy model and a more re- 
alistic and physically-motivated model for reionization. The 
first one, the tanh model, is expressed in terms of two pa- 
rameters, namely the two main characteristics of the over- 
all reionization process, its midpoint z r and span Az; the 
particular form of the model is merely mathematically con- 
venient, with no real physical significance, and it restricts 
reionization to be smooth and symmetric about its mid- 
point. The second model, the CDM model, is based on the 
standard understanding of galaxy formation within CDM- 
dominated halos. It assumes a fixed overall ionizing effi- 
ciency Mon (number of ionizing photons per baryon), a den- 
sity clumping factor C and a minimum halo circular velocity 
V c for galactic halos, and it yields reionization models with 
up to 3 parameters (C is relatively minor and we typically 
held it fixed in the fitting) . Unlike the tanh model, the CDM 
model is asymmetric, with the exponentially increasing halo 
abundance leading to an acceleration of reionization in its 
later stages. 

Despite the fact that the tanh model is a simple parame- 
terization that has often been used in the literature, we have 
shown that it leads to substantial systematic errors if it is 
assumed when fitting a 21-cm signal that is described by the 
CDM model. In particular, the best-fit z r in the tanh model 
is underestimated typically by 5-10%, while Az is underes- 
timated by tens of percent. However, a sufficiently sensitive 
experiment (e.g., with an integration time ii n t > 20 hours 
for the case of a foreground polynomial of degree N po \ y — 3) 
would be able to discriminate between the CDM and tanh 
models based on the x 2 value of the best-fit model. 

Our main result is a detailed plot of the detection limits 
of global 21-cm experiments (Figure [12}. We find that the 
CDM model can produce quick reionization scenarios (with 
a redshift span Az ~ 1) if feedback makes large halos dom- 
inate, which then requires a high ionizing efficiency in these 
halos (see also Figure Some of these realistically possible 
models can be ruled out with 50-hour global 21-cm exper- 
iments even in the pessimistic case where a polynomial of 
degree JV po iy = 12 is required for removing the foregrounds 
(or other systematic effects). If somewhat more ambitions 
experiments are achievable, then a broad range of scenarios 
up to z r ~ 12 can be probed within our CDM model. The 
smooth and symmetric tanh model is more difficult to dif- 
ferentiate from the foreground polynomial, and it requires 
greater integration times and lower -/V po i y in order to rule 
out for similar reionization characteristics. 

Our conclusions are generally optimistic in terms of the 
possibility for global 21-cm experiments to reconstructing 
the reionization history and constrain the properties of the 



ionizing sources. In particular, one-year EDGES observa- 
tions may allow a remarkably precise reconstruction. How- 
ever, the polynomial degree N po \ y that is required for remov- 
ing the foreground and systematic effects plays an important 
role. In the most optimistic case, where N po i y — 3 suffices 
over the entire frequency range of 100-250 MHz, 1% errors 
on z r are achievable; they require t- m t — 51 hours with the 
tanh model or 29 hours with the CDM model (in each case 
with Az = 2 as defined in that model). These, of course, 
are only statistical errors, while we have shown that there 
can be much larger systematic errors if the assumed reion- 
ization model cannot reproduce the real 21-cm signal from 
reionization. 

Indeed, our results merit some caution, since our in- 
vestigation indicates that a broader range of flexible and 
realistic models of reionization should be studied before we 
can be confident that the results are robust. For instance, 
the parameters of the model (JVi on , V c , and C) could change 
with redshift due to evolving feedbacks such as metal en- 
richment or the effect of photoheating on suppressing gas 
accretion onto galaxies in the reionized regions (this effect 
is large if reionization is initially dominated by relatively 
small hal os). Such an ev olution could, e.g., be parameter- 
ized as in lBarkanal (120091 '); while such a complication of the 
model would no doubt lead to serious partial degeneracies 
among the parameters, hopefully the main characteristics 
of reionization would remain measurable at high accuracy. 
Another possible complication that could be added is the in- 
creasing effect of recombinations near the end of reionization 
due to the optical depth of dense clumps w ithin the then- 
large HII bubbles (jFurlanetto fe Ohl [20051 ). We note that 
some of these effects, such as the feedback and the increas- 
ing recombinations, should slow the progress of reionization 
as it nears its end, thus rounding the shape of the CDM 
model and perhaps reaching a result that is more similar 
to the tanh model. We plan to study these more realistic 
models. 

We note that in the radiometer equation we neglected 
the effects of the instrumental response (or bandpass) on 
both the foreground and the cosmological signal (i.e., we as- 
sumed a flat bandpass filter). Indeed, in a real observation 
this would raise three issues, but it should be possible to deal 
with them (Judd Bowman, personal communication). First, 
our theoretical model must be convolved with the instru- 
mental response, in order to compare to the observations. 
This can be done if the bandpass is known with sufficient 
accuracy (~1% in current experiments and expected to im- 
prove). The foregrounds must also be convolved with the 
response, but as long the instrumental response is smooth, 
i.e., does not introduce sharp frequency features, the fore- 
grounds can still be removed by low-order polynomial fit- 
ting. Finally, another effect of a less than perfect response is 
a decrease of the sensitivity of the observation with respect 
to a flat receiver response. We can simply compensate for 
this effect by slightly increasing the observing time (by a 
few tens of percent for current experiments). 

Recently, iBowman fc Rogers! (|2010T ) reported substan- 
tial new results from their upgraded EDGES experiment. 
Taking only the cleanest data out of their observational run, 
they had a thermal noise level equivalent to tint = 0.8 hour 
(with our idealized assumptions) within their 100-200 MHz 
band. They also found it necessary to use iV po i y = 5 in or- 
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der to remove the foreground and systematics, and reach the 
thermal noise level, within 20 MHz sub-bands in their spec- 
trum. This is still somewhat worse than even our conserva- 
tive iVp i y = 12 case ove r the full 100-250 MHz ran ge. Even 
with these limitations, iBowman fc Rogers! l|2010T ) reached 
an observational milestone, namely the first direct obser- 
vational limit on the rapidity of cosmic reionization. In par- 
ticular, using the tanh model, they set a 95% confidence 
lower limit of Az > 0.06 for the duration of the reionization 
epoch. Global 21-cm experiments are still in their infancy 
and clearly have a quite promising future. 
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